First-principles characterization of ferromagnetism in N-doped SrTiC>3 and BaTiOa 
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The spin-polarization and magnetic coupling character of N-doped SrTiC>3 and BaTiC>3 are studied 
through first-principles calculations. The substitutional N doping at O sites leads to a half-metallic 
property and produces a magnetic moment of 1.0 /^s. The magnetic interaction between the nearest 
and next-nearest N dopants results in a strong ferromagnetic coupling. When the distance between 
the N dopants is larger than 7 A, the ground state of the system tends to be paramagnetic. A 
nitrogen-concentration threshold to produce the ferromagnetism is estimated. The calculated results 
give a good explanation for the experimentally observed ferromagnetism in N-doped SrTiC>3 and 
BaTi0 3 . 



Diluted magnetic oxides have attracted much atten- 
tion owning to their fundamental physics and promis- 
ing applications in the spintronics.[l| With respect to 
the transition-metal doped oxides, an unexpected room- 
temperature ferromagnetism (RTFM) has recently been 
observed in C-doped ZnO.Q Soon after, lots of first- 
principles theoretical calculations have been done to pre- 
dict the possible ferromagnetism in C (N)-doped oxides, 
such as ZnCsd Q Ti0 2 @-[3 and Sn0 2 ,[l and the 
RTFM in these C (NVdoped oxides was confirmed by the 
later experiments. |Rj-[l4l| Recently, the RTFM was found 
in N-doped SrTi0 3 (STO)[l| and BaTi0 3 (BTO). OH It 
is well-known that the perovskite STO and BTO are im- 
portant ferroelectric materials, [l?) and thus it is of great 
interest to realize the coexistence of ferroelectricity and 
ferromagnetism in the two oxides for their multiferroic 
applications. [l8| Liu et al.'s experimental studies on the 
N-doped STO show that the ferromagnetism only ap- 
pears in the sample with a relatively high doping con- 
centration, and N doping leads to an insulator-metal 
conducting phase transition. (l5| As a result, it is spec- 
ulated that there might exist a nitrogen-concentration 
threshold, only above which can the N-doped STO and 
BTO exhibit the long-range ferromagnetism. To address 
this question and understand the origin of its ferromag- 
netism, we studied the spin-polarization and magnetic 
coupling character of N-doped STO and BTO on the 
basis of the first-principles calculations. The theoreti- 
cal results show that the N doping leads to the typi- 
cal half-metallic properties, and produces a strong fer- 
romagnetic coupling through the magnetic orbital inter- 
action between the nearest and next-nearest N dopants. 
A nitrogen-concentration threshold to produce the long- 
range ferromagnetism is also estimated. 

Various doped structures with N atoms at O sites are 
modeled using the 3x3x3 supercell containing 135 atoms 
based on the cubic perovskite STO and BTO, as shown 
Fig. 1. The projector augmented wave (PAW) poten- 
tials are used for electron-ion interactions and general- 




' kcsong.yang@gmail.com 



FIG. 1. (Color online) 3x3x3 supercell employed to simulate 
N-doped STO and BTO. The O atoms in pink color labeled 
by 0-9 are the sites to be replaced with N atoms. The x-, 
y- and z- axes are along the crystallographic a-, b- and c- 
directions, respectively. 



ized gradient approximation parameterized by Perdew- 
Burkc-Ernzcrhof (PBE) is used for exchange-correlation 
functional. [ll| The cut-off energy of 400 eV for the plane- 
wave basis set and a 2 x 2 x 2 k-point grid centered at 
T point are used. [20[ The convergence threshold for self- 
consistent iteration is set at 10 -6 eV, and all the atomic 
positions are fully optimized until all components of the 
residual forces are smaller than 0.01 eV/A. All the cal- 
culations are performed using the Vienna ab-inito simu- 
lation package (VASP).[2l|, 

To understand the origin of the spin-polarization in 
N-dopcd STO and BTO, we firstly studied one- N- atom 
doped model, which is constructed by replacing one O 
atom at the site using one N atom in the supercell (sec 
Fig. 1). The calculated total density of states (TDOS) 
and partial DOS (PDOS) for N-doped BTO and STO 
are shown in Fig. 2a (a') and 2b (b'), respectively. It 
shows that the substitution of N for O introduces some 
spin-polarized impurity states in the band gap, and they 
mainly originate from N 2p orbitals. The Fermi level 
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FIG. 2. (Color online) The total DOS and partial DOS plot 
for one-N-doped (a) STO and (b) BTO. The dotted line in- 
dicates the Fermi level at eV. 



is pinned in the middle of these band-gap states, indi- 
cating that N-doped STO and BTO exhibit the typi- 
cal half-metallic character. This is good agreement with 
the previous theoretical calculation |23[ and Liu et al.'s 
experiment that N-doped STO shows the conducting 
property. (l5| With respect to the local Cartesian coor- 
dinate defined in Fig. 1, the up-spin and down-spin N 
2p x orbitals are occupied, and the up-spin N 2p y and 2p z 
orbitals are also occupied. However, the Fermi level is 
pinned in the middle of the down-spin N 2p y and 2p z or- 
bitals, i.e., half of the down-spin N 2p y and 2p z orbitals 
are occupied. This indicates that one N atom at an O 
site generates a magnetic moment of 1.0 /is, and the N 
dopant exists as a N 2_ (s 2 p 5 ) anion. In addition, non- 
spin-polarizcd calculations for one-N-atom doped STO 
and BTO show that the spin-polarized state is more sta- 
ble than the non-spin-polarized state by about 124 and 99 
meV, respectively, thus indicating that the ground state 
of N-doped STO and BTO is magnetic. 

To examine the type and strength of magnetic coupling 
between the two N dopants, nine inequivalent structural 
configurations are modeled using the 135-atom 3x3x3 
supercell, which are constructed by replacing two O 
atoms using two N atoms at positions (0, 1), (0, 2), (0, 
3), (0, 4), (0, 5), (0, 6), (0, 7), (0, 8) and (0, 9), re- 
spectively, as shown in Fig. 1. For convenience, we use 
(i, j) to denote the structure in which the O atoms at 
the positions (i, j) are replaced by N atoms. To examine 
the stable magnetic coupling type, both the ferromag- 
netic (FM) and the antiferromagnetic (AFM) coupling 
between the moments of two N dopants are considered 
for each structure. In TABLE I and II, we summarize the 
calculated results for N-doped STO and BTO. For each 
structure, the optimized N- • -N distance and the rela- 
tive energy AE, the magnetic energy E mag =E FM -E AFM , 
the total magnetic moment (M) under the ferromagnetic 
alignment and the calculated Curie temperature (Tc) are 
listed. For the (i, j) configuration, the relative energy AE 
is defined as the energy difference between the energy of 
the lower-energy (either FM or AFM) and the total en- 
ergy of the (0, 1) structure, i.e., AE=E(i, ])-E(0, 1). 
The Tc is calculated based on the mean-field theory and 



TABLE I. Calculated results for N-doped STO: the opti- 
mized N- • -N distance (djv-Jv) (A), the relative energies AE 
(eV), magnetic energy E mag =E F M-EAFM) (meV), the total 
magnetic moment (M) (/is /cell) under ferromagnetic align- 
ment and Curie temperature (Tc) (K) calculated for the (i,j) 
structure of the two-N-atom doped STO. 
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TABLE II. Calculated results for N-doped BTO. 
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Heisenberg model, i.e., k B T c =(2/3)E mag .\M, [25[ Here 
Emag is the magnetic energy. 

Table I and II both show that the relative energies AE 
of various (i, j) configurations are nearly uniform, indi- 
cating that the N dopants in STO and BTO do not have 
the energetic preference to form a cluster. For the (0, 
1), (0, 2) and (0, 3) structures, a strong FM coupling 
occurs between the magnetic moments of two N dopants. 
This is consistent with the ex per imentally observed ferro- 
magnetism in N-doped STOpl] and BTO.pH Further- 
more, the highest calculated Tc of N-doped STO and 
BTO both exceed 1000K, much higher than that of C- 
doped ZnO (about 400 K).@, Q Therefore, N-doping is 
an effective route to obtain high-Tc half-metallic ferro- 
magnetism in STO and BTO. For (0, 5) structure, a 
stronger FM coupling takes place in N-doped STO than 
in doped BTO. As the N- ■ -N distance increases, i.e., 
for the (0, 4), (0, 6), (0, 7), (0, 8) and (0, 9) struc- 
tures, their magnetic energies decrease a lot, and even 
become zero when the N- • -N distance exceeds about 
7A. This indicates that N-doped STO and BTO tend 
to be paramagnetic when the distance between the two 
N dopants is long enough. It means that an effective 
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(a) (b) 

FIG. 3. (Color online) Calculated spin density distribution of 
two-N-atom doped STO for (a) (0, 1) configuration and (b) 
(0, 4) configuration. 

magnetic coupling interaction only occurs between the 
two nearest and next-nearest neighbor N dopants. As 
a result, we speculate that there exists a percolation 
threshold of nitrogen-concentration. [26| and if the mag- 
netic interaction range is limited to the next-nearest N 
dopants, a minimal nitrogen-concentration of about 11.1 
% is essential to produce the long-range ferromagnetism. 
This gives a good explanation for the experimental fact 
that the ferromagnetic hysteresis loop only appears in N- 
doped STO with a relatively high concentration (see Fig. 
4 in Ref.[H). 

To understand the ferromagnetic coupling mechanism 
in N-doped STO and BTO, the (0, 1) and (0, 4) struc- 
tures of N-doped STO are chosen as examples to show 
their three-dimensional spin density distributions, which 
are presented in Fig. 3. The spin density is mostly con- 
tributed by the N dopants as well as their nearest O 
atoms. For (0, 1) structure, the magnetic orbitals of the 
two N dopants have a substantial direct overlap, while 
for (0, 2), (0, 3) and (0, 5) structures, though not shown 
here, an indirect magnetic orbitals overlap through their 
common O atoms is found. Therefore, a direct or indirect 
magnetic orbitals overlap between the two N sites could 



be responsible for the FM coupling. [f| With respect to 
the other remaining configurations, the N- • -N distance 
in (0, 4) structure is much shorter, and to clearly indi- 
cate the magnetic coupling interaction between two N 
dopants, its spin density distribution is shown in Fig. 
3b. In this configuration, the two N dopants induce the 
spin-polarization of the nearest O atoms in the different 
plane (Nl is in xz plane, and N2 in yz plane). Hence, the 
magnetic orbitals of the two N dopants could only have a 
weak overlap through only one common O atom (labeled 
as 01 in Fig. 3b), thus resulting in a weak magnetic 
interaction. For the other structures, the long N- • -N dis- 
tances determine that the magnetic orbital overlap be- 
tween the two N dopants is small or even zero, and thus 
their magnetic ground states tend to be paramagnetic. 

In summary, our first-principles calculations for N- 
doped STO and BTO indicate that the N doping in- 
troduces spin-polarized impurity states in the band gap, 
generating a magnetic moment of 1.0 /is, and leads 
to the half-metallic property. The substantial mag- 
netic interaction between the nearest and next-nearest 
N dopants leads to the strong ferromagnetic coupling, 
but the ground state tends to be paramagnetic when the 
distance between the N dopants is larger than about 7 
A. A minimal nitrogen-concentration to produce long- 
range ferromagnetism is estimated. These results give a 
good explanation for the experimental phenomena, and 
provide a possible route to realize their multiferroic ap- 
plications. 
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